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£NJ , Abstract 

£C) • The total cross section for the production of massive quarks in electron positron 

annihilation can be predicted in perturbative QCD. After expansion in m 2 /s the 

■ quartic terms, i.e. those proportional to to 4 /s 2 , are calculated up to order for 
vector and axial current induced rates. Predictions relevant for charm, bottom and 
top quarks production are presented. The corrections are shown to be comparable 
to terms of order a s and a 2 . As a consequence, the predictions exhibit a sizeable 
dependence on the renormalization scale. The stability of the prediction is improved 
and, at the same time, the relative size of the large order terms decreases by replacing 

■ the running mass fh(ji) with the scheme independent invariant one m. By combining 
CLi these results with the prediction for massless case and the quadratic mass terms the 

■ cross section for massive quark production at electron positron colliders is put under 
D ' control in order from the high energy region down to fairly low energies. 

& : 

■ The total cross section for hadron production in electron-positron annihilation is one of 
the most fundamental observables in particle physics. For energies sufficiently far above 
threshold it can be predicted by perturbative QCD, and it is well accessible experimentally 
from threshold up to the highest energies of LEP and a future linear collider. It allows 
for a precise determination of the strong coupling a s and, once precision measurements at 
different energies are available, for a test of its evolution dictated by the renormalization 
group equation. In many cases the cms energy is far larger than the quark masses which 
motivated the original calculations to be performed in the idealized case with the masses 
set to zero from the start. In this limit, the results of O (o^) Q and O (a^) [|| were 
obtained more than two and nearly one decade ago, respectively (for a review see ||). 

However, in a number of interesting cases quark masses do play an important role Q. For 
Z decays into bottom quarks the mass effects have to be included as a consequence of the 
extremely precise measurements. Bottom quark production at lower energies is affected 
from production threshold up to a few tens of GeV. Other cases of interest Q[| are charm 
production between roughly 5 and 10 GeV and, last but not least, top quark production 
at a future linear collider [pi. 
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In two-loop approximation the full mass dependence of the cross section has been evaluated 
since long and, exploiting the optical theorem, both real and imaginary parts are avail- 
able 0. In three loop approximation (O (a 2 )) the corresponding results were obtained 
only recently. Two methods have been used for this purpose: the first one is based on 
the evaluation of a large number of terms for the Taylor series of the polarization function 
H(q 2 ) at q 2 = and an appropriately chosen analytic continuation ||]. The second one is 
based on the application of the large momentum expansion (see || and references therein), 
which provides an expansion of H(q 2 ) in powers (m 2 /q 2 ) n , modulo logarithms |[[C|| . From 
the comparison between full result and expansion one learns that the first few terms of 
the high energy expansion provide a remarkably good description of the full result, from 
high energies down to values of 2m/y / s = 0.65 - 0.75. The existence of the four quark 
threshold at = 4m and of a corresponding branching point for the polarization func- 
tion suggests that the high energy expansion diverges for ^fs below 4m. Nevertheless, 
numerical studies |L(]] as well as qualitative arguments demonstrate that the sum of the 
first two or three terms of the expansion can be trusted even down to 3m or even, with 
some optimism, 2.5m (where m stands for the pole mass). 

These considerations pave the way to a prediction of R{s) including the quark mass de- 
pendence to O (a^) along the following route: In addition to the massless result the m 2 /s 
terms of O (af) have been calculated for the absorptive part nearly a decade ago 11]. They 



were obtained by reconstructing the logarithmic a^m 2 /s terms of H(q 2 ) from the full three 
loop O (a 2 m 2 /s) result of |[^] with the help of the renormalization group equations. These 
are sufficient to calculate the m 2 /s terms of the imaginary part in the time-like region. 
The result is cast into a particularly compact form once expressed in terms of the running 
mass m(fi 2 ), with the 't Hooft scale \i set to y/s throughout fi"3|| , since this choice elimi- 
nates all terms oc ln(s/m 2 ). A generalization of this approach has been formulated for the 



quartic mass terms in [14,15] and was originally adopted for the calculation of a 2 m 4 /,s 2 
terms. It is based on the operator product expansion and the usage of the renormalization 
group equation to again construct the logarithmic terms of the polarization function. In 
addition to the anomalous mass dimension and the /^-function the anomalous dimensions 
of the operators of dimension four are required in appropriate order. 

This method also allows to determine the a^m 4 /s 2 terms. The calculation can be reduced 
to the evaluation of massless propagators and massive tadpole integrals, both at most in 
three loop approximation. Details of the calculation will be given elsewhere [O]. Com- 
pared to the case of the m 2 / s terms an important difference arises: Even adopting the MS 
definition (at scale /i 2 = s) for the quark mass, logarithmic terms oc ln(s/m 2 ) remain in 
order a 2 and above. The power of these logarithms, however, is reduced by one. In fact, 



these logs may be also summed up (see Refs. [17, 1 



To predict the R ratio, it is well justified to consider only one quark as massive and to 
neglect the masses of the lighter quarks. The heavier quarks decouple (apart from the 
tiny axial-vector singlet contribution — see below) and can therefore be omitted in the 
calculations. The massive quark will generically be denoted by Q in the following, while 
q refers to all the lighter quarks. 

In the energy region where quark mass effects are relevant, charm and bottom production 
is solely induced by the electromagnetic vector current. Top quark production, however, 
which sets in above 350 GeV, receives additional contributions from the axial current. Both 
vector and axial-vector will be considered separately in the following. The prediction for 
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both of them is conveniently split up as follows: 

R (v) = 



E2 , 2 («) i (v) 
v q r q + v rW +r, v 



sing 



(1) 



and similarly for the axial- vector part (v — ► a). The sum runs over all massless quark 
flavors. v q /Q {a q /q) is the coupling constant to the light/heavy quarks in the vector (axial- 
vector) case. For low energies only the electromagnetic current contributes, v — * e and 
a — > 0. For high energies, both the electromagnetic and neutral current pieces are relevant 
and have to be included with appropriate weights (see, e.g., fill ). r q and Tq represent 
the non-singlet contributions arising from diagrams where the external currents are linked 
by a common quark line. They originate from two different types of diagrams: For r q the 
external current couples to massless quarks; the massive quark then only appears through 
its coupling to virtual gluons. r q is the same for external vector and axial-vector currents. 
On the other hand, r^q ^ corresponds to diagrams where the external current couples to 



the massive quark. 
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, finally, comprises massless and massive singlet contributions, 



where either of the external currents is coupled to a separate closed quark line. 



Both r q and r. 



(v/a) 



are written as series in rrig/s, where s is the cms energy and rriQ is the 
MS mass of the quark Q: 



r q = ro+ r qj2 + r qA + 



J v / a ) — fr\ —I— r^l^ 4- r^l^ + 



2) 



Q,2 <2,4 

ro denotes the massless approximation, while r q ^ n and are the mass terms of order 

rriq . If not stated otherwise, the renormalization scale fi 2 = s is adopted below. 

Denoting rif the number of active flavors, the massless terms are given by (all the following 
formulas are valid up to O (af ) unless indicated otherwise) 
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The quadratic mass corrections, separated according to @), read ||IT, 18|: 
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The quadratic mass terms in r q contribute in O (o? s ) and higher only [11|. Finally, we 
present the quartic mass corrections up to O (a 3 ): 
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The singlet contributions in the vector numerically small: 
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They do not receive m 2 corrections [11|. 

The singlet contributions in the axial-vector case are exceptional in the sense that for 
bottom quark production the top quark does not decouple, i.e., the contributions where 
the top quark couples to one of the external currents are not suppressed by powers of 



1/m 2 . The corresponding corrections up to O (m 2 «g) have been computed in |20|-23|. 
They turn out to be small, so we will neglect them in the following. 

Considering the axial vector singlet contribution to top quark production, on the other 
hand, one should take into account the full top-bottom doublet, because the axial anomaly 
cancels in this combination. However, this strategy induces completely massless final 
states without any top quarks. At order a 2 the expansion for the contributions from the 
full top-bottom doublet starts with the to 6 terms ||. The separate piece from massless 
cuts is known in analytical form [20 and has to be subtracted if one wants to obtain the 
contribution from the final states with top quarks. The complete result is given in Ref. ||. 
Thus the singlet contribution to tt production is given in expanded form by 
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where the weak coupling a 2 = a 



-a t a b 



1/4 has been pulled out for clarity. The 



separate contribution from the massless cuts is not yet available in order a^. In principle 
it should be subtracted from the following result in order to arrive at the production rate 
for top quarks: 

/ o\ 2 



sing 



a. 



TT 



mr 



380 
— C3 + 



1520 



/40 lou 

Cs + n f [ — Cs - — Cs 



160 
9" 



6 





5. 


6. 
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8. 
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10. 


10.5 


«i 4j ( S ) 


0.212 


0.2 


0.192 
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0.174 
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0.829 


0.803 


0.784 


0.769 


0.756 


0.745 


0.74 





11. 12. 14. 16. 20. 




0.174 0.171 0.165 0.161 0.153 
3.61 3.57 3.5 3.44 3.35 





420. 460. 500. 




0.097 0.0961 0.0952 
154. 153. 152. 



Table 1: Running coupling and masses at different scales. Matching for a s is 
performed at the values m(m) given in (|14"|). 



(373.116 - 13.0918n/) . (13) 

The prime indicates that this expression still contains the contributions from purely mass- 
less final states, as mentioned before. We explicitly denote the mass by m t here in order to 
recall that, on the one hand, this result applies only to top production, and on the other 
hand, the full (t, b) doublet has been taken into account. 

Let us now investigate the numerical significance of mass effects for the charm and bottom 
case. If not stated otherwise, we will adopt the following input data: 

M z = 91.19 GeV, a s {M 2 z ) = 0.118 , 

m c {m c ) = 1.2 GeV , m b {m b ) = 4.2 GeV , m t (m t ) = 165 GeV . (14) 

The running of the quark masses and the coupling constant to the scale -^/i is performed 
with three loop accuracy. Since we are working in the MS scheme, we have to take into 
account the matching conditions for as when going from rif to n/ ± 1. We perform the 
matching from nj = 5 to nj = 4 at mb(mb) and from n/ = 5 to rif = 6 at mt{mt). For 
illustration, some values for a s (s), m c (y / s), mb(y / i), and mt(y / i) are displayed in Table [j]. 
The charm, bottom, and top masses are defined for rif = 4, 5, and 6, respectively. As 
stated above, the approximation is expected to become unreliable in the threshold region, 
that means below 5.5 — 6 GeV for charm, 12 — 12.5 GeV for bottom, and 420 — 450 GeV for 
top contributions, which justifies the choice for the cms energies in the figures presented 
below. 

In Figs. [l]-|| we display separately the contributions for the massless case (Fig. |]) and for 
the m 2 and m 4 terms, including successively higher orders in a s . While Fig. || shows the 
effects of the diagrams with light quarks coupling to the external current, Figs. || and |] 
correspond to the case where the massive quark couples directly to the external current. 
The axial contribution is presented for the top quark only. 

The variation of the prediction with the renormalization scale \x is shown in Figs. [| to |8| 
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Table 2: r^> for yfi = 10.5 GeV. 



As discussed already in earlier papers, the size of the higher order corrections decreases 
quickly with increasing order in a s , as far as the massless approximation is concerned 
(Fig. [H) . This is reflected in the stability of the result under variation of the renormalization 
scale [i by a factor between 1/2 and 2 as displayed in Fig. ||[ The m 2 terms in r q (Fig. ||| (a), 
(c), (e)) arise for the first time in order af. As a consequence this prediction exhibits a 
strong n dependence as displayed in Fig. ^ (a), (c), and (e). However, this term is typically 
around 10 -4 and thus irrelevant for all practical purposes. 

The quartic terms in r q which are displayed in Figs. || (b), (d), and (f) contribute in O (a 2 ) 
and higher. Terms of order m 4 a^ and m 4 a 2 are of comparable magnitude, not only for the 
low energy case but even at the highest energy, whence the prediction for the m 4 terms in 
r q must be considered as uncertain within a factor two. This is reflected in the strong /i 
dependence of the result shown in Fig. ||| (b), (d), and (e). 

No reduction of this variation is observed when moving from the ai to the a^ calculation. 
Nevertheless, the prediction for the total cross section is not seriously affected by this 
instability since these terms are again of order 10 -4 only and thus far below the foreseeable 
experimental precision. 

The dominant quartic mass terms are obviously expected from tq, the part where the 
massive quark is coupled to the external current. For comparison we also discuss the 
quadratic terms which are known since long and are shown in Fig. ^. The vector cur- 
rent induced piece is displayed in Fig. || for c, b, and t quarks, the axial piece is shown 
for t quarks only. Terms of increasing order in a s decrease in magnitude and apparent 
convergence is observed. This welcome behavior is reflected by the improved stability of 
higher orders under variations of \x (Fig. |7|) . The behavior of the quartic terms which are 
the theme of this paper is more problematic (Fig. Q). The bulk of the result is given by 
the Born approximation. The order a s -, a 2 -, and (^-corrections are significantly smaller 
than the leading terms. However, with increasing order they do not decrease but remain 
roughly comparable in magnitude. This is reflected in the strong fi dependence displayed 
in Fig. |8|. Depending on the choice of ji, the relative size of the three corrections varies 
drastically. 

Nevertheless, the higher orders are small compared to the m 4 Born terms and, even more 
important, their instability does not affect the stability of the overall prediction for R 
resulting from the sum of the different terms. Accepting as an estimate of the uncertainty 
of rQ t 4 its variation with fj, between yfs/2 and 2-y/s, r Cj 4 at 6 GeV varies by ±0.0005, 
at 14 GeV by ±0.0016. For r^4 the variation is negligible. 

These considerations demonstrate that a prediction for R has been obtained which is valid 
in order a$ and includes mass terms in an expansion up to order m 4 . 

From a pragmatic point of view the smallness of the m 4 terms, in particular of their QCD 
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corrections, allows to ignore their apparent instability. Nevertheless, one may also attempt 
to arrive at a more stable result for the quartic terms by replacing the running mass by 
the invariant one (fh) through 

m(a) = rh exp f da ^"! . , (15) 
J ap{a) 

where a = a s (fi 2 )/iT. (3(a) and 7 m (a) are the renormalization group functions governing 
the running of the strong coupling constant and the quark mass: 

^ 2 -^ja = a(3(a) , /i 2 -^jm = mj m (a) . (16) 
d/i z d// z 

Their perturbative expansion is known up to O (a*) both for /3(a) and 7 m (a) |2q| . 
The integral in (jl^) is solved perturbatively, and the resulting expression for R is re- 
expanded in a s up to the third power. For the scale invariant mass we assume the values 
rh c = 2.8 GeV, rhb = 15.3 GeV, and fat = 1076 GeV which corresponds to solving (|i~5| ) 
w.r.t. fh for fi = m at three-loop order and using the numerical values of Eq. (]i~4|). The 
result for the \i dependence of Tq \ in this approach is shown in Fig. ^. 

The variation of the 0(a^) prediction with /i is reduced. For example, for = 6 GeV 
the O(a^) prediction in Fig. 8 varies between -0.10 • 10" 2 and -0.20 • 10~ 2 with a central 
value of —0.18 • 10~ 2 at fi = yfs, compared to a range from —0.13 • 10~ 2 to —0.20 • 10 -2 
with a central value of —0.19 • 10~ 2 in Fig. 9. At the same time one observes a reduction 
of the higher order terms compared to the Born approximation and the 0(a s ) prediction. 
For fi = yfs the results within the two approaches are quite close, which gives additional 
confidence in the reliability of these numbers. 



Summary 

The total cross section for the production of massive quarks in electron positron annihila- 
tion can be predicted in perturbative QCD. After expansion in m? / s the quartic terms, i.e. 
those proportional to m 4 /s 2 , were calculated up to order for vector and axial current 
induced rates. Predictions relevant for charm, bottom and top quark production were 
presented. The a? corrections were shown to be comparable to terms of order a s and a 2 . 
As a consequence, the predictions exhibit a sizeable dependence on the renormalization 
scale. Adopting instead of the MS scheme a framework where the running mass fh(ij) is 
replaced by the invariant mass fh, the stability of the prediction is improved and, at the 
same time, the relative size of the large order terms decreases. Obviously, an improved 
understanding of the origin of these large corrections would be highly desirable. 

Combining these results with the massless prediction and the quadratic mass terms we 
have demonstrated that the cross section for massive quark production at electron positron 
colliders is under control in order from the high energy region down to fairly low 
energies. 



Remark 



The results of this paper are available in Mathematica format at 

http : //www-ttp .physik.uni-karlsruhe . de/Progdata/ttp00/ttp00-09/. 
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Figure 1: The massless approximation tq in three different energy ranges, relevant 
for (a) charm, (b) bottom and (c) top production. 
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Figure 2: Mass corrections to the non-singlet contribution of r q for c, b, and 
t production (1. ,2., and 3. row, respectively), arising from diagrams where the 
external current couples to massless quarks only. Left column: quadratic, starting 
to be non-zero in O(a^); right column: quartic mass corrections, starting to be 
non-zero in O {a 2 s ). 
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Figure 3: Quadratic mass corrections (oc m?) to the non-singlet contribution of r<g 
for Q = c,b, t, arising from diagrams where the external current couples directly to 
the massive quark. Upper row: vector current for c and b quarks; lower row: vector 
and axial currents for t quark. 
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Figure 4: Quartic mass corrections (oc to 4 ) to the non-singlet contribution of rq 
for Q = c,b, t, arising from diagrams where the external current couples directly to 
the massive quark. Upper row: vector currents for c and b quarks; lower row: vector 
and axial currents for t quark. 
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Figure 5: Variation of tq with \x. s is fixed to values relevant for the production of 
(a) charm, (b) bottom, and (c) top quarks. The abscissa ranges from [i = \fsj2 to 
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Figure 6: Variation of r q p (left column) and r 9; 4 (right column) with fi. The 
first, second, and third row correspond to charm, bottom, and top production, 
respectively. 
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Figure 7: Variation of ^ with \i. First row: (a) charm, (b) bottom production 
for vector case; second row: (c) vector and (d) axial-vector case for top production. 
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Figure 8: Variation of with \i. Conventions as in Fig. [?]. 
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Figure 9: Variation of Tq 4 a with using the invariant mass rh. The conventions 
are the same as in Fig. |7| and |8[ 
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